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ABSTRACT 

A characterization of the optical turbulence vertical distribution (CTjy profiles) and 
all the main integrated astroclimatic parameters derived from the C^- and the wind 
speed profiles above the site of the Large Binocular Telescope (Mt. Graham, Arizona, 
US) is presented. The statistic includes measurements related to 43 nights done with 
a Generalized Scidar (GS) used in standard configuration with a vertical resolution 
AH~1 km on the whole 20 km and with the new technique (HVR-GS) in the first 
kilometer. The latter achieves a resolution AH^20-30 m in this region of the atmo- 
sphere. Measurements done in different periods of the year permit us to provide a 
seasonal variation analysis of the Cfj . A discretized distribution of Cjj useful for the 
Ground Layer Adaptive Optics (GLAO) simulations is provided and a specific analysis 
for the LBT Laser Guide Star system ARGOS (running in GLAO configuration) case 
is done including the calculation of the 'gray zones' for J, H and K bands. Mt. Gra- 
ham confirms to be an excellent site with median values of the seeing without dome 
contribution e = 0.72", the isoplanatic angle 8q = 2.5" and the wavefront coherence 
time To = 4.8 msec. We find that the optical turbulence vertical distribution decreases 
in a much sharper way than what has been believed so far in proximity of the ground 
above astronomical sites. We find that 50% of the whole turbulence develops in the 
first 80±15 m from the ground. We finally prove that the error in the normalization 
of the scintillation that has been recently put in evidence in the principle of the GS 
technique, affects these measurements with an absolutely negligible quantity (0.04"). 
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1 INTRODUCTION 

The Mt. Graham International Observatory (MGIO) is lo- 
cated on Mt. Graham (32°42'05" N, 109°53'31" W), Arizona 
(US) and hosts three telescopes: the Vatican Advanced Tech- 
nological Telescope (VATT - D = 1.83 m) the Heinrich Hertz 
Submillimiter Telescope (SMT - D = 10 m) and the Large 
Binocular Telescope (LBT - two D = 8.4 m dishes that, when 
working in an interferometric configuration, can achieve the 
resolution of a telescope with a 23 m pupil size). The study 
and characterization of the optical turbulence (OT) distribu- 
tion in space and time is fundamental for the ground-based 
astronomy in the visible up to the near-infrared range to 
design adaptive optics systems and to optimize their perfor- 
mances. The vertical distribution of the OT (i.e. the C% 
profiles) is the parameter from which all the integrated as- 
troclimatic parameters derive from. 

In this paper we present a study based on measurements 



of the C% profiles related to 43 nights and obtained with 
a Generalized Scidar (GS) placed at the focus of the VATT 
on the Mt. Graham summit, around 250 m from the LBT. 
Figure 4 in Egner & Masciadri (2007) shows the relative 
position of the two telescopes. It is worth to note that the 
primary mirror of the LBT is located ~ 35 m above the 
dome of the VATT. 

The scientific motivations for such a long-term site test- 
ing campaign are: 

(1) To collect an as rich as possible statistical sample of 
optical turbulence (OT) vertical distribution (C% profiles) 
to be compared with simulations obtained with the atmo- 
spherical mesoscale model Meso-Nh with the aim to validate 
the model above Mt. Graham. This is a key milestone for 
the ForOT projectQ whose final goal is to predict the optical 
turbulence above astronomical sites (Masciadri 2006). The 
measurements from a vertical profiler such as a Generalized 
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Scidar are crucial for the validation of such a kind of mod- 
els. It is our interest to collect a heterogeneous sample of 
measurements taken in different periods of the year and dif- 
ferent turbulence conditions in a way to better control the 
model behaviors under different conditions and to better val- 
idate the model itself. The atmospherical models has been 
used for the first time to reconstruct and characterize the 
C% profiles by Masciadri et al. (1999a, 1999b). Since then 
many progresses have been achieved by our group: the model 
has been applied to different astronomical sites in a simple 
monomodel configuration (Masciadri et al. 2001, Masciadri 
& Garfias 2001) and more recently in a grid-nesting configu- 
ration (Lascaux et al. 2009) , a new calibration technique has 
been proposed (Masciadri & Jabouille 2001) and statistically 
validated (Masciadri et al. 2004) and the first application of 
the Meso-Nh model as a tool of turbulence characterization 
has been presented (Masciadri & Egner 2006). However, so 
far we always could access to GS measurements concentrated 
in a well defined period of time. This series of site testing 
campaigns with a GS at Mt. Graham aimed to comply with 
the necessity to diversify the measurements sample. 

(2) To provide a characterization of all the most impor- 
tant integrated astroclimatic parameters above the site of 
the Large Binocular Telescope (LBT) and verify if any ev- 
ident changes are observed with respect to previous re- 
sults obtained on a three times smaller sample (Egner et 
al. 2007). This study can be, therefore, important to con- 
firm/reject/refine those conclusions. We remind that the GS 
is a manually operated instrument and the statistical rich- 
ness of GS measurements can not be compared to that of 
monitors such as the DIMM and the MASS that are rou- 
tinely run above many observatories. However, these mon- 
itors can provide only integrated values (as is the case for 
the DIMM) or turbulence distribution with very low verti- 
cal resolution AH ~ h/2 (as it is the case for the MASS). 
The GS remains therefore a unique instrument for a detailed 
analysis of the vertical distribution of the OT in the whole 
troposphere. 

(3) To provide an as rich as possible statistical sample of the 
high resolution vertical distribution (AH = 200-250 m and 
25-30 m) of the optical turbulence in the first hundreds of 
meters up to 1 km to support the feasibility studies of new 
generation instruments for the LBT, such as the LBT Laser 
Guide Stars system ARGOS (Rabien et al. 2008), that, in its 
first baseline, is planned to work with a GLAQj configura- 
tion. The GLAO efficiency, indeed, strongly depends on the 
turbulence distribution and strength near the ground. The 
study we intend to perform can be achieved thanks to a 
new technique that has been recently proposed called High- 
Vertical Resolution Generalized Scidar (HVR-GS - Egner & 
Masciadri 2007) that aims to reconstruct the optical turbu- 
lence vertical distribution in the first kilometer above the 
ground with a resolution 4 up to 10 times higher than what 
has been done so far with standard vertical profilers such as 
the GS (typically AH ~ 1km). In Egner & Masciadri (2007) 
the validity of this technique has been proved and it is now 
our intention to characterize the turbulence distribution in 
statistical terms, to verify how the turbulence decreases in 
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the first hundreds of meters above the ground and provide 
inputs to test the ARGOS performances. 

In Section [5] we briefly review the principle of the GS 
and HVR-GS techniques. In Section [3] we present the site 
testing campaign data-set that has been analyzed and pre- 
sented in a preliminary form in Stoesz et al. (2008). In Sec- 
tion|3]we present an exhaustive statistical analysis of all the 
most important integrated astroclimatic parameters derived 
from the C% profiles and characterizing Mt. Graham. In 
Section [S]we quantitatively discuss the turbulence distribu- 
tion of the OT in the whole troposphere in application to the 
Adaptive Optics. We provide a composite profiles distribu- 
tion on the whole troposphere with the vertical resolution 
required for the characterization of the LBT-LGS system 
running in GLAO configuration (ARGOS) i.e. (AH ~ 200 
m in the first kilometer and AH ~ 1 km above 1 km). In 
Section [6] we characterize the C% profiles at standard and 
high-vertical resolution and their seasonal variation. Section 
[7]summarizes the conclusions. In Annex we briefly show that 
the error in the normalization of the scintillation of the GS 
technique, put in evidence by Johnston et al. (2002) and 
Avila & Cuevas (2009), produces effects absolutely negligi- 
ble on these measurements. From a general point of view, 
we will show that GS measurements obtained with a pupil 
size D > 1.5 m and a binary separation 9 < 8 arcsec, are 
affected by this error for less than a few hundredths of an 
arcsec. 



2 INSTRUMENTS: GS AND HVR-GS 

Two instruments have been used for this study: the GS and 
the HVR-GS. We used the GS as developed by McKenna 
et al. (2003). The Scidar technique has been originally pro- 
posed by Rocca et al. (1974) and Vernin & Azouit (1983) and 
relies on the analysis of the scintillation images generated by 
a binary in the pupil plane of a telescope. The Scidar tech- 
nique (called Classic Scidar) is insensitive to the turbulence 
near the ground. More recently Fuchs et al. (1998) proposed 
a generalized version of the Scidar (called Generalized Sci- 
dar) in which the detector is virtually conjugated below the 
ground permitting to extend the measurements range to the 
whole atmosphere (~ 20-25 km). The GS principle has been 
later put in practice by several authors above different astro- 
nomical sites (Avila et al. 1997, 2004, Klueckers et al. 1998, 
McKenna et al. 2003, Fuensalida et al. 2004, Egner et al. 
2007, Garcia-Lorenzo et al. 2009). The GS is based on the 
observation of binaries with a typical separation 8 within 
(3-10) arcsec, the binary magnitude mi,ni2 < 5 mag and 
A(mi : 2) < 1 mag. When two plane wavefronts propagating 
from a binary meet a turbulence layer located at a height 
h from the ground, they produce on the detector plan, op- 
tically placed below the ground at about h 9S = 3 km, two 
scintillation maps made by a set of characteristic shadows 
appearing in couple separated by a distance d that is ge- 
ometrically related to the position of the turbulence layer 
as d = 6- (h + h 9S ). The calculation of the auto-correlation 
(AC) of the scintillation map produces the so called 'triplet'. 
The central peak is located in the centre of the AC frame, 
the lateral peaks are located at a symmetric distance d from 
the centre. The amplitude of the later peaks is proportional 
to the strength of the turbulence of the layer located at the 
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height h weighted by the scintillation that such a layer pro- 
duces on the detector. 

In a multilayer atmosphere different turbulent layers 
produce triplets with lateral peaks located at a different dis- 
tance d* from the centre of the AC frame. To monitor the 
whole troposphere the pupil size needs to be large enough 
(D > 1.5 m) to avoid that one of the shadows that form a 
couple (or both of them) falls outside the pupil of the tele- 
scope. In the AC frame the triplets are all placed along the 
direction of conjunction of the binary. To retrieve the C% 
profile, the central peak of the triplet in which some sources 
of noise such as the photon noise are present, is firstly elim- 
inated. Finally the C% is obtained inverting the Fredholm 
equation that, in the LBT-GS, is done using the conjugated 
gradient (Egner et al. (2007)) algorithm. The vertical res- 
olution AH(h) of the GS technique depends on our ability 
in discriminating two different later peaks and it is propor- 
tional to the Fresnel Zone (FZ) (Vernin & Azouit, 1983): 

= 0.78 -y/y^ = 0_78_FZ 
9 9 

For typical values of the observable binaries with a standard 
GS the typical vertical resolution is AJf(O) ~ 1 km. 

The HVR-GS technique has been introduced recently 
(Egner & Masciadri, 2007) and aims to measure the C% 
profiles with a high vertical resolution (AH ~ 25-30 m) in 
the first kilometer above the ground. We briefly summarize 
here the main concepts on which this technique is based on. 
If we abandon the idea to monitor the whole 20 km and 
we limit our attention to the first kilometer we can easily 
increase the vertical resolution of a GS up to a factor 4 
using the standard GS technique based on the calculation 
of the autocorrelation (AC) obtained from the scintillation 
maps of binary stars having a separation 9 around four times 
larger than the typical separation used for the standard GS 
technique. This is not enough, however, to achieve resolu- 
tion of the order of 25-30 m because we are fundamentally 
limited by the Fresnel Zone size. However, if we use simulta- 
neously the autocorrelation (AC) and the cross-correlation 
maps (CC) taken with 20-40 msec time lag, the triplets in 
the CC frames are not aligned anymore on the same direc- 
tion identified by the binary separation but the central peak 
of each triplet is located on the direction of the wind speed 
of each turbulent layer distributed in the troposphere and 
the distance from the centre of the CC frame is proportional 
to the wind speed of the same turbulent layer. Under this 
assumption it is possible to prove (Egner & Masciadri 2007) 
that the vertical resolution Ah p i x is finally limited by the 
accuracy with which we can estimate the position of the lat- 
eral peaks in the CC frames that is smaller than the pixel 
size projected on the pupil (A x p i x ) and it is equal to: 

Ah pix = 0.56^P (2) 

9 

For 9 = 30 arcsec and A x p i x = 7 mm we retrieve a 
typical vertical resolution A h P i X ~ 25-30 m. In conclu- 
sion, the main HVR-GS concept consists in: (1) taking a 
wide binary of the order of 30-35 arcsec monitoring the first 
kilometer and (2) treating simultaneously the AC and CC 
frames. From a practical point of view, the C% retrieved 
from the AC frames is characterized by an energy that is 
redistributed in a set of thinner layers within the first kilo- 



Obscrving Runs Nights 



27 April 2005 1 

19-24 May 2005 6 

6-15 December 2005 5 

27 May 2007 - 3 June 2007 8 

16-28 October 2007 13 

23 February 2008 - 3 March 2008 10 



Table 1. Observing runs at Mt. Graham. 

meter whose vertical resolution is of the order of 25 m. The 
AC frames are fundamental to re-normalize the total energy 
in the first kilometer. The final high vertical resolution C% 
profile is retrieved from: 

h max 

j C%,Ac{h) ■ dh = f 'scale ' ^ ' C 2 N , CC (hi) (3) 

-Ah max (0)/2 1 

where h max is the height of the highest detected layer, the in- 
ferior limit of the integral Ah max (0) /2 takes into account the 
vertical resolution at h=0. f SC aie is a factor that takes into 
account the decorrelation of the central peak of the triplets 
having a V > ± AV (AV = 0.2-0.8 m-s" 1 ) with respect 
to the zero velocity triplet. It has been measured (Egner & 
Masciadri, 2007) that i S caie in the range AH = 1 km can 
be considered the same for all the thin turbulent layers. We 
note that, using the HVR-GS we have to consider a rate of 
rejection of frames because of the wind fluctuations that can 
produce a spreading of the central peak of the triplets. The 
rate of rejection can be more or less conservative, depending 
on the constraints that the user wishes to introduce in this 
analysis. Anyway, in Egner & Masciadri (2007) it has been 
proved that no bias in the measurements is introduced if we 
reduce the statistical sample. This can be explained with the 
fact that the conditions that facilitate the spreading of the 
central peaks are not necessarily correlated to typical bad 
or good seeing. We refer the reader to that paper for further 
details. We highlight that two slightly different approaches 
have been proposed recently in the literature to increase the 
vertical resolution near the ground: the HVR-GS (Egner & 
Masciadri 2007) and LOLAS (Avila et al. 2008, Avila & 
Chun, 2004). The main difference from the point of view of 
the principle is the following: the HVR-GS uses a known 
instrument (the Generalized Scidar) but with a new tech- 
nique. LOLAS uses a new instrument but with the same 
technique of the Generalized Scidar. The readers can find 
in the corresponding papers the details (and technicalities) 
related to both. 



3 SITE TESTING CAMPAIGNS: DATA-SET 
STATISTIC 

We collected and analyzed so far observations related to 43 
nights (Table [TJ. As already noted the GS is manually oper- 
ated therefore the site testing campaigns have been sched- 
uled respecting the periods in which the VATT was shut 
down (July and August) and trying to cover as many dif- 
ferent periods of the year. Table [2] shows the code used to 
identify the typology of the data-set: 'GS' indicates the stan- 
dard GS measurements extended on ~ 20 km (43 nights), 
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Sample 


Nights 


Measurements 


Hours 


Resolution 


'GS' 


43 


16657 


163 


AH(0) ~ 1km 


'WB' 


15 


3659 


6.2 


AH(0) ~ 200 m 


'HVR-GS' 


15 


2812 


5.1 


AH(0) ~ 25 m 



Table 2. Classification of the GS campaign measurements. 'GS': 
C%, profiles retrieved from the standard GS. 'WB': Ci. profiles 



retrieved from the AC frames of the GS measurements obtained 
with wide-binaries. 'HVR-GS': Cfj profiles retrieved from the 
wide-binaries auto-correlation (AC) and cross-correlation (CC) 
frames following the technique described in Egner & Masciadri 
(2007). 



Name 


"J2000 


<5j2000 


mi 


iri2 


e 








(mag) 


(mag) 


(arcsec) 


7 Ari 


01 54 


+ 19 17 


4.5 


4.6 


7.6 


Castor 


07 35 


+ 31 53 


1.9 


3.0 


4.4 


7 Leo 


10 20 


+ 19 50 


2.4 


3.6 


4.7 


it Boo 


14 41 


+ 16 25 


4.9 


5.8 


5.5 


S Ser 


15 34 


+ 10 32 


4.2 


5.2 


1.1 


95 Her 


18 02 


+ 21 36 


4.9 


5.2 


6.5 


P Cyg 


19 31 


+ 27 58 


3.2 


4.7 


35.3 



Table 3. Binary stars observed with the GS and HVR-GS. 



'WB' indicates the C% retrieved from the autocorrelation 
(AC) frames associated to wide binaries (AH ~ 200-250 m) 
(15 nights) and 'HVR-GS' the C% obtained with high ver- 
tical resolution (AH ~ 25-30 m) (15 nights) in which both 
the autocorrelation and cross-correlation frames have been 
treated. The samples of the three categories are differently 
rich because we started to use the HVR-GS more recently 
and this new technique has a higher rate of measurements re- 
jection. We note that, to avoid biases in the estimates in the 
HVR-GS, it has been decided to discard from the statistic 
all doubtful cases characterized by the presence of clouds or 
cirrus. A method that we called 'normalization' (and it will 
be described later) has been applied to the sample of high 
resolution measurements (15 nights). Among other advan- 
tages it permits us to provide a turbulence budget represen- 
tative of the whole sample of 43 nights. Basically, the mor- 
phology of the turbulence energy distribution (shape of the 
C% versus the height) is retrieved from the observation of 
wide-binaries for 15 nights and with the 'normalization' pro- 
cedure the turbulence energy of the first kilometer detected 
with the standard GS is redistributed in thin turbulent lay- 
ers according to the profile morphology reconstructed with 
the CC frames. The selected binary stars for the standard 
GS and HVR-GS techniques are reported in Table They 
are substantially the same indicated in Egner et al. (2007) 
and Egner & Masciadri (2007) but we eliminated from the 
sample those binaries with magnitude larger than 6 mag (i.e. 
118 Tau has been discarded) to avoid potential biases due 
to a too weak intensity. 



4 INTEGRATED ASTROCLIMATIC 
PARAMETERS 

The seeing, the isoplanatic angle 9o, the wavefront coherence 
time to and the equivalent velocity Vo are defined as: 



ro = 



0.423 



C 2 N {h)dh 



-3/5 



e = 0.98— 

ro 



% = 0.057 ■ A B/5 



h s/3 C%(h)dh 



-3/5 



V = 



JV(h) 5/3 C 2 N (h)dh 
o 

oo 

/ C 2 N {h)dh 



3/5 



(4) 
(5) 
(6) 

(7) 
(8) 



r = 0.31 — 

Vo 

Figures [1] shows the cumulative distribution of the as- 
troclimatic parameters (seeing, seeing in the free atmosphere 
calculated for h > 1 km, isoplanatic angle and wavefront 
coherence time) calculated for the total 43 nights and for 
the (April-September) and (October-March) periods that we 
will simply call summer and wintei0. Table[j]summarizes the 
median (50 th ), first (25 th ) and third (75 th ) quartiles for the 
three main integrated astroclimatic parameters calculated 
for the following groups: the whole sample of 43 nights, the 
summer and the winter time. 

A composite wind speed profile has been used to calcu- 
late the median wavefront coherence time ro. Below 2 km the 
wind speed retrieved from the GS has been used; above 2 km 
the wind speed profile retrieved from the European Centre 
for Medium-Range Weather Forecast (ECMWF) analyses 
extracted in the nearest grid point (32.75°N, 110.00° wjj to 
the Mt. Graham (~ 11.5 km northwest of the summit). Due 
to the fact that the wind speed vertical profiles retrieved 
from the ECMWF analyses are calculated at the synoptic 
hours (00:00, 06:00, 12:00, 18:00) a temporal interpolation 
of adjacent synoptic hours wind speed have been performed 
to better represent the wind speed in the local temporal 
range in which C% measurements have been done. It has 



3 As shown in Table[T] measurements do not cover all the months 
of the year with an exact uniform distribution. For example there 
are no measurements in the July-September period. However we 
are interested in putting in evidence macroscopic differences be- 
tween the two extreme seasons and we observed that the mor- 
phology of the average Cfj profile and the turbulence strength 
was substantially the same including or not the data related to 
the October month. We therefore decided for this division instead 
of a more reductive (April- June) and (December-March) division 
that would have imply to study a less statistically representative 
sample. 

4 It is worth to highlight that the analyses from the ECMWF are 
available, at present, with a horizontal resolution of 0.25 degrees. 
Calculations done in Egner et al. (2007) have been obtained with 
data extracted at latitude (33° N) and a horizontal resolution of 
0.5 degrees. 



Optical turbulence vertical distribution at Mt. Graham 5 



Median Seeing Median Free Atmosphe 




(") (msec) 



Figure 1. Cumulative distributions of four integrated astroclimatic parameters. Top- left: the total seeing (including the dome contri- 
bution). Top-right: the seeing in the free atmosphere (h > 1 km). Bottom-left: the isoplanatic angle (0q). Bottom-right: the wavefront 
coherence time (to). Thick lines: the whole sample. Dotted lines: summer time. Thin lines: winter time. 





Total 






Summer 






Winter 






Parameter 


25 th 


50 th 


75 th 


25*'' 


50 th 


75 th 


25*'' 


50 th 


75 th 


e (arcsec) 


0.65 


0.95 


1.34 


0.53 


0.61 


0.72 


0.89 


1.19 


1.50 


9q (arcsec) 


1.6 


2.5 


3.6 


3.1 


3.8 


4.5 


1.4 


2.0 


2.7 


to (msec) 


2.7 


4.8 


8.7 


6.4 


10.1 


14.6 


2.5 


3.8 


6.2 


V (ms" 1 ) 


5.1 


7.2 


9.3 


3.7 


5.1 


7.8 


5.9 


7.7 


9.6 



Table 4. Median, first and third quartiles values of the main integrated astroclimatic parameters above Mt. Graham (43 nights): seeing 
in the total troposphere, isoplanatic angle, wavefront coherence time, integrated equivalent wind speed. 



been already shown (Egner et al. (2007)) that this is the 
best method to treat the wind speed to calculate a reliable 
tq. The ECMWF analyses do not well reconstruct the oro- 
graphic effects produced on the atmospheric flow at the top 
of the summit. Moreover, the ECMWF grid points spaced 
by 0.25 degrees correspond to locations with lower altitudes 
than the summit altitude Ho and, as a consequence, the 
ECMWF wind speed calculated at an height equal to Ho 
is, in general, larger than the wind speed measured on the 
summit. On the other side, measurements of the wind speed 
with GS implies a great number of rejected frames and it 
is frequently difficult to retrieve a profile extended on the 
whole troposphere. The composition method for the wind 
speed revealed, therefore, to be the best solution for the 
calculation of to (Egner et al. 2007). 

Looking at Fig[T]and Table [4] a clear seasonal variation 
appears evident for all the integrated astroclimatic parame- 
ters. In May 2005 and May 2007, in both years, 10-15 days 
of extremely good seeing and large isoplanatic angle features 
has been observed. Such a long time of extremely good con- 
ditions in the same period of the year indicate that this is, 
highly probably, among the best periods of the year for tur- 



bulence conditions at Mt. Graham. This result is coherent 
with the typical weather conditions at synoptic scale in this 
region and in this season that are characterized by weaker 
wind speed in the high troposphere. A low probability to 
trigger optical turbulence in the high atmosphere confirmed 
by a weaker C% strength (see Section [6}, is coherent with a 
large median 9q and a small median seeing in the free atmo- 
sphere (Fig. [TJ in the summer seasons. If we also take into 
account the typical weaker equivalent wind speed (Eq(4]| in 
summer with respect to winter time (see Tableland Section 
[6]| we can explain the typical larger value of the median to 
in this season (see Eq(8]). 

Finally, to quantify the contribution of the seeing 
provided only by the atmosphere, the dome seeing (ej), 
calculated with the method described in Egner et al. (2007) 
and Avila et al. (2001), has been subtracted from the 
total seeing (e). The method consists in discriminating 
the triplets located at the ground (H=H(0)±AH) with 
a velocity V = ± AV from those located at the same 
height and having a velocity different from zero (V > AV). 
The velocity resolution AV of our system is 0.8 to 0.2 
m/s per pixel with time lags within the range (10 to 40 
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Median dome seeing 




0.5 



(") 



Figure 2. Cumulative distribution of the dome seeing for all of 
the 43 nights (thick line), the summer (dotted line) and winter 
(thin line) period. 



turbulence conditions, the pitch size of the adaptive optics 
system. A more detailed discussion on that will be done later 
but it is, of course, evident that it can be very useful for the 
GLAO applications (1) to know the budget of turbulence 
energy developed in the high part of the atmosphere i.e. how 
much turbulence remains not corrected by a GLAO system 
and (2) to better estimate the size of the gray zone and the 
vertical distribution of the turbulence inside the gray zone. 
The gain of a GLAO system depends, indeed, mainly on 
these two issues. From the standard GS measurements it is 
trivial to calculate the cumulative distribution of efa for 
h > 1 km (Figure [TJ. The same calculation is done for the 
whole year, the summer and winter time. We calculated a 
median value sfa = 0.39 arcsec for h > 1 km, (ef^,,,™ = 
0.31 arcsec, SFA,win = 0.44 arcsec). For those cases in which 
Umax = 1 km the efa represents the portion of turbulence 
that is not corrected at all. In the next session we will deal 
about the portion of turbulence developed above X m, where 
X < 1km extremely important for GLAO applications. 



Parameter 


Egner et al. (2007) 


Masciadri et al. (2009) 


e included dome 


0.80 


0.95 


e without dome 


0.68 


0.72 


00 


2.71 


2.5 


TO 


3.6 


1.8 



Table 5. Comparison of values obtained in the study of Egner 
et al. (2007 ) - 16 nights and Masciadri et al. (2009) - 43 nights. 
Seeing (with and without dome contribution) and £>n are in arcsec, 
the tq in msec. 



ms). Figure [2] shows the cumulative distribution of the 
dome seeing calculated for the whole sample and the two 
seasons. The median value of the 'dome seeing' is e a = 
0.52 arcsec (Fig(5J. For the first time it has been observed 
an interesting seasonal trend in the dome contribution 
that certainly deserves a carefully investigation in the 
future. Knowing that the median seeing in the whole 
atmosphere (included the dome seeing) is £=0.95 
arcsec and that e a = 0.52 arcsec, it follows that 
the median seeing related to the whole atmosphere 
without the dome contribution for the richest statis- 
tic we collected so far (43 nights) is £tot= 0.72 arcsec. 

Table [S] reports the difference in the median estimates of 
the principal integrated astroclimatic parameters obtained 
above Mt. Graham with samples that are differently rich by 
a factor ~ 3. No major differences can be highlighted with 
exception of the tq that appears larger of about 1 msec with 
respect to the paper of Egner et al. (2007). 

To complete the analysis of the integrated astroclimatic 
parameters it is worth to remind that, in the context of 
GLAO systems (Tokovinin, 2004), the atmosphere can be 
divided in three vertical slabs: a region near the ground (h 
< H min ) in which the turbulence is totally corrected, an 
intermedium region (H m i„ < h < Y{ max ), called gray zone 
in which the turbulence is partially corrected and a region 
that covers the rest of the troposphere (h > H maa; ) in which 
the turbulence is not corrected any more. The values of H m i„ 
and H max depend on the wavelength, the field of view, the 



5 COMPOSITE PROFILES FOR ADAPTIVE 
OPTICS APPLICATIONS 

The 'composite profiles ' are used to represent efficiently and 
in statistical terms the vertical distribution of the optical 
turbulence in a discretized number of layers particularly 
suitable for AO simulations. The method is commonly used 
by many authors in the field (Tokovinin & Travouillon 2006, 
Egner et al. 2007, Stoesz et al. 2008, Masciadri et al. 2009, 
Chun et al. 2009) and it consists in identifying a finite num- 
ber of vertical slabs covering the whole troposphere (~ 20 
km) with their correspondent value of seeing (Ae^ or AL) 
so that the total turbulence integrated on the troposphere 
is conserved. J is defined as: 

oo 

J = J C%{h)dh (9) 



and it is related to the seeing as: 
J = 9-KT n -A 1/3 -e 5/3 (10) 

where A is expressed in meters, e in arcsec and J in m 1 ' 3 . 
The turbulence in the free atmosphere (h > 1 km) and in the 
boundary layer (h < 1 km) is treated in an independent way 
to permit to study different combinations of probabilities for 
the OT vertical distribution. Table |S] reports the AJ; values 
in the range h > 1 km calculated at different heights and 
obtained from the C% associated to the rn related to the 
20-30 % of its cumulative distribution ('good' case), to the 
45-55 % ('typical' case) and to the 70-80 % ('bad' case). 
Measurements from the sample 'GS' (see Table [5J are used. 

The C% profiles retrieved from the 'WB' sample char- 
acterized by a Ah ~ 200 m near the ground (H < 1 km) have 
a suitable vertical resolution to calculate the composite pro- 
files in this vertical range for applications to ARGOS (field 
of view 6 — 4 arcmin) because the turbulence developed be- 
low H mln = AX/(2-6>) ~ 200 m (AX = 0.5 m is the pitch size 
i.e. the projection of the actuator of the deformable mirror 
on the pupil of the telescope) is resolved by the instrument. 
Basically we do not need a higher vertical resolution for this 
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application. We therefore first calculated a similar distribu- 
tion of Table [6] for h < 1 km and obtain a temporary table 
that we call table T similar to Table [6] As already antici- 
pated in Stoesz et al. (2008) and Masciadri et al. (2009), to 
take into account the different statistic obtained with the 
standard 'GS' and wide binary sample 'WB' and to take 
into account the quantitative information of the turbulence 
present in the first kilometer provided by the whole sample 
of 43 nights, each number of this temporary table T has to 
be multiplied by the correction factor f 9 ;: 



reported in Table [7| to finally obtain the composite dis- 
tribution for h < 1 km (Table [SJ). The numerator of f 9 ; is 
the seeing measured in the first kilometer from the ground 
from the 'GS' sample. The denominator is the seeing mea- 
sured in the first kilometer from the ground from the 'WB' 
sample. Table [8] reports the composite profiles for h < 1 
km equivalent to 43 nights. It includes the dome contribu- 
tion frequently preferable for AO simulations. In the last 
row the J values in the case in which the dome contribution 
is subtracted are reporte The multiplication by the f 9 ; 
factor (method that we call normalization) offers the great 
advantage to retrieve the spatial distribution of the turbu- 
lence in the first kilometer using the HVR-GS technique (or 
the 'WB' as is this case) and to use the quantitative tur- 
bulence energetic budget of data extracted from the stan- 
dard GS technique that can be done on a richer statistical 
sample (43 nights). We overcome therefore the intrinsic lim- 
itation of a high number of rejected frames typical of the 
cross-correlation technique. This method also is completely 
insensitive to any bias potentially introduced by the normal- 
ization of the autocovariance of the scintillation maps of the 
GS put in evidence by Johnston et al. (2002) and Avila & 
Cuevas (2009) (see Annex B for demonstration). The median 
value of the composite profiles (Table [6] and Table [8] central 
columns) permits us to calculate the percentage of the tur- 
bulence developed above different heights h with respect to 
the turbulence developed in the whole troposphere (Table[9]- 
left side) as well as the percentage of turbulence developed 
in the (0, h) range with respect to the turbulence devel- 
oped in the whole troposphere (Table fright side) . Table 
1101 shows the same calculation obtained in case the median 
dome seeing is subtracted. 

The composite distribution of the C% permits also to 
quantify the extent of the gray zone. We note, firstly, that 
Umax =t/0 (where r is the residual wavefront coherence size 
after GLAO correction i.e. the residual FWHM) depends on 
the wavelength. Using the C% composite profiles extended 
on 20 km as an input (Tableland Table |SJ), with the com- 
bination 'good'-'good', 'typical'- 'typical' and 'bad'-'bad', in 
the boundary layer and in the free atmosphere, the ARGOS 
simulations provide the residual values r (L. Busoni, private 
communication). From these values of r that we reported in 
Table [JTJ we can retrieve the values of H maa; for the wave- 
lengths in the near-infrared range: J, H and K (Table ITT)) . 

5 We note that the median seeing of Table [4] is not exactly the 
same as the seeing retrieved by the composite profiles because 
the former treats a set of individual seeing values that are pro- 
portional to J 3 / 5 . 



Bins 


'Good' 


'Typical' 


'Bad' 


(m) 


J (m 1 / 3 ) 


J (m 1 / 3 ) 


J (m 1 / 3 ) 


14000-20000 


7.61e-15 


1.02e-14 


1.91c-14 


12000-14000 


5.10c-15 


9.65e-15 


1.57e-14 


10000-12000 


7.24e-15 


1.27e-14 


2.89c-14 


8000-10000 


8.52e-15 


1.72e-14 


3.22e-14 


6000-8000 


6.45e-15 


1.05e-14 


1.77c- 14 


4000-6000 


9.50c-15 


1.23e-14 


2.10c-14 


3000-4000 


9.18c-15 


1.14e-14 


1.45e-14 


2000-3000 


1.97c-14 


3.39e-14 


4.02c-14 


1500-2000 


6.98c-15 


1.52e-14 


2.86c-14 


1000-1500 


5.28e-15 


1.50e-14 


2.86e-14 



Table 6. Composite profiles for h > 1 km. In the first column 
the boundaries of the vertical slabs. In each column is reported 
the value of J that is proportional to the integral of the optical 
turbulence in the correspondent vertical slab. These composite 
profiles are statistically representative for 43 nights. 



Ground Layer Seeing 


'Good' 


'Typical' 


'Bad' 


GS 


0.55 


0.81 


1.16 


WB 


0.62 


0.75 


0.88 




0.82 


1.1 


1.6 



Table 7. f g ; factor for the 'Good', 'Typical' and 'Bad' distribu- 
tion. 



Knowing that H m i„ = 200 m, as explained previously, we 
deduce that the 'gray zone' extends in the (200 m - 378 
m) range and assume its smallest value when we observe 
in J band and we consider the 'bad'-'bad' case (75% case). 
It extends in the (200 m - 3777 m) range and assumes its 
largest value when we observe in K band and we consider 
the 'good'-'good' case (25% case). Fig. [3] shows the Umax 
(blue, red and yellow lines) for different field of view in the 
case the residual FWHM of GLAO simulations has been ob- 
tained with the central column of Table [6] and Table [8] i.e 
(case 50%). 



Bins 


'Good' 


'Typical' 


'Bad' 


(m) 


J (m 1 / 3 ) 


J (mV3) 


J (m 1 / 3 ) 


900-1000 


4.35e-15 


7.06e-15 


1.33e-14 


800-900 


2.48e-15 


3.53e-15 


6.29e-15 


700-800 


3.30e-15 


6.04e-15 


1.04e-14 


600-700 


5.72e-15 


1.09e-14 


2.05e-14 


500-600 


4.30e-15 


7.70e-15 


1.38e-14 


400-500 


4.21e-15 


9.60e-15 


1.86e-14 


300-400 


2.05e-14 


4.26e-14 


8.26e-14 


200-300 


6.15e-15 


1.14e-14 


2.53e-14 


100-200 


2.03e-14 


4.14e-14 


8.11e-14 


0-100 


1.79e-13 


3.27c-13 


6.17e-13 


0-100 


3.4e-14 


8.5e-14 


2.20e-13 



Table 8. Composite profiles for h < 1 km after the 'normaliza- 
tion' for the f g i factor. These composite profiles are statistically 
representative for 43 nights. In the last line are reported the J 
values obtained without the dome contribution (median values: 
e d,25 = 0-35 arcsec, £^50 = 0.52 arcsec, £^,75 = 0.70 arcsec). 
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75% 


FWHM 


H-max 


50% 


FWHM 


Umax 


25% 


FWHM 


Umax 




(arcsec) 


(m) 




(arcsec) 


(m) 




(arcsec) 


(m) 


J 


0.43 


378 


J 


0.30 


567 


J 


0.18 


945 


H 


0.37 


606 


H 


0.20 


1123 


H 


0.11 


2042 


K 


0.25 


1208 


K 


0.13 


2324 


K 


0.08 


3777 



Table 11. Values of K ma x calculated for 9 = 4 arcmin and different residual FWHM. The FWHM values are obtained with the GLAO 
simulations using as inputs Table [6] and Table [8] H m i n = AX/26 ~ 200 m for all the wavelengths. 



h 


J[h,top}/Jtot 


h 


J[0,h]/Jt°t 


(m) 


(%) 


(m) 


(96) 


1000 


25 


1000 


75 


800 


26 


800 


74 


600 


29 


600 


71 


400 


31 


400 


69 


300 


38 


300 


62 


200 


40 


200 


60 


100 


47 


100 


53 



D = 8 m, pitch = 0.5m, f wh m_K = 0. 1 3", fwhm_H = 0.2", fwhm_J=0.3" 



Table 9. Dome seeing included - Left: Percentage of turbu- 
lence developed above the height h with respect to the turbulence 
developed on the whole troposphere. Right: Percentage of turbu- 
lence developed between the ground and the height h with respect 
to the turbulence developed in the whole turbulence. Second and 
fourth columns are obviously complementary. 



h 


J[h,top]/Jtot 


h 


J[0,fi] /Jtot 


(m) 


(%) 


(m) 


(%) 


1000 


40 


1000 


60 


800 


43 


800 


57 


600 


47 


600 


53 


400 


52 


400 


48 


300 


63 


300 


37 


200 


66 


200 


34 


100 


77 


100 


23 



Table 10. Dome seeing subtracted - Left: Percentage of tur- 
bulence developed above the height h with respect to the turbu- 
lence developed on the whole troposphere. Right: Percentage of 
turbulence developed between the ground and the height h with 
respect to the turbulence developed in the whole turbulence. Sec- 
ond and fourth columns are obviously complementary. 



6 OPTICAL TURBULENCE VERTICAL 
DISTRIBUTION: C% 

6.1 GS: vertical distribution on the whole 
troposphere 

Figure [4] shows the median C% profile obtained with the 
whole data-set of 43 nights, the summer and the winter pe- 
riods. The morphology of the vertical distribution of the op- 
tical turbulence (C N profile) shows that the greatest tur- 
bulence contribution develops in the first kilometer above 
the ground. Between 1 and 10 km we observe a set of mi- 
nor peaks changing their position and strength during the 
year. At around 10 km we observe the typical secondary 
C% peak developed at the jet-stream level. For what con- 
cerns the seasonal variation we observe that the ground layer 
bump, responsible for most of the turbulence budget, shows 




Figure 3. Extent of the 'gray zone' i.e. H m ; n < h < Umax for 
different fields of view and wavelengths in the case of median 
distribution (50%). H ma:E (colored lines) is calculated for J (blue), 
H (red) and K (yellow) band. H m i n (black line) is the same for all 
the wavelengths. The pupil size is D = 8 m and the pitch size AX 
= 0.5 m. The FWHM is equivalent to r i.e. the residual wavefront 
coherence size after correction. 



a clear seasonal trend indicating larger turbulence strength 
in winter than in the summer period. In the free atmosphere 
we observe the interesting effect of the secondary C% peak 
located at 10 km in winter time that shifts to higher heights 
(~ 14 km) and is characterized by a weaker strength in sum- 
mer time. The latter effect (that we call 'a effect') has been 
put in evidence the first time by (Masciadri & Egner 2006) 
above a different site (San Pedro Martir) with simulations 
provided by a mesoscale atmospherical model. Above Mt. 
Graham the jet-stream C% peak appears to have a similar 
shift of ~ 4 km toward higher heights in summer and the 
C% peak is located at roughly the same absolute height 
from the ground with respect to San Pedro Martir (Masci- 
adri & Egner, 2006). At that time we had no measurements 
extended on different periods of the years above San Pedro 
Martir to retrieve a seasonal trend. Avila et al. (2004) re- 
ferred just to the spring and therefore they could not put 
in evidence any seasonal trend. It is worth to note that this 
was, to our knowledge, the first time that an atmospherical 
model could put in evidence new insights before measure- 
ments showing that simulations can be a valuable tool to 
investigate the nature of the turbulence by itself. 

More recently (Els et al. 2009) MASS measurements 
done above San Pedro Martir and extended on yearly time 
scale have been published as part of the TMT site testing 
project (Schoeck et al. 2009). The MASS vertical resolution 
is much lower than the GS one but some information, useful 
in this context, can be retrieved. The layer at 8 km (the 
nearest to the jet-stream level) is unfortunately not perfectly 
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centred on the jet-stream and considering a resolution of 4 
km at this height the absolute values of the peak-to-peak 
seasonal variation can be smoothed out. However, if we do 
not take care about the absolute estimate of the amplitude 
peak-to-peak variation for which the MASS is not the most 
suitable instrument, we can say, in any case, that the MASS 
8 km layer can provide a qualitative seasonal trend near the 
jet-stream level. In Els et al. (2009) - Fig. 5 we can observe 
above San Pedro Martir a seasonal effect of the turbulence 
strength similar to what has been obtained above the same 
site with simulations performed with atmospheric models 
(Masciadri & Egner 2006) and to what has been observed 
above Mt. Graham i.e. an increase of the turbulence strength 
in the local winter at the jet-stream level. A similar trend 
is evident in Els et al. (2009)-Fig.5 also for all the other 
astronomical sites. It is worth to note that in the sites in 
the south hemisphere, the seasonal trend of the 8 km layer 
is inverted i.e. we found the maximum values in proximity of 
the local winter. This effect seems therefore confirmed so far 
above all the astronomical sites for which measurements are 
available. However, the MASS is not the suitable instrument 
to put in evidence the more complex effect ('a effect') of 
the vertical shift of the jet-stream C% peak towards higher 
heights in summer as is resolved by the GS. Also the lack 
of seasonal trend of the layer located at 16 km in the Els 
et al. (2009) paper does not mean that there is no seasonal 
variations at this height but simply that this instrument can 
not resolve seasonal variations since its vertical resolution is 
~ 8 km at h = 16 km. We conclude therefore that results 
obtained at San Pedro Martir and at Mt. Graham with a 
GS are an evidence that the 'a effect' exists and it is not 
typical of a specific site. 

Similar behaviors of the C% observed above different 
sites is promising for describing a 'universal physical model' 
able to explain what is the origin of the a effect. Masciadri 
& Egner (2006) proposed an explanation. 

Starting from the assumption that the development of 
the optical turbulence depends on the gradient of the wind 
speed and the potential temperature a parallel analysis of 
these two elements all along the troposphere in the two sea- 
sons can provide insights on the seasonal variations trends. 
Above San Pedro Martir no great difference have been iden- 
tified in the potential temperature and its gradient between 
the summer and winter. On the other hand a substantial 
decrease of the gradient of the wind speed and its strength 
has been identified at a jet-stream level in summer. This 
fact could explain the decrease of the turbulence strength in 
summer at a jet-stream level. Besides, the rapid inversion of 
the wind speed gradient at around 15-16 km could explain 
the production of the optical turbulence at these heights 
that appears as a shift of the jet-stream C% peak towards 
higher heights. 

Figure [5] shows the median potential temperature and 
wind speed profiles above Mt. Graham calculated in the two 
periods (winter and summer) in which C% measurements 
have been collected. The potential temperature is retrieved 
from the ECMWF analyses in the same grid point we dis- 
cussed in Section [4] The wind speed profiles are obtained 
with the composite procedure described in Section [4] i.e. 
GS measurements have been considered for h < 2 km and 
ECMWF analyses for h > 2 km. We observe a not negli- 
gible wind speed increase at the jet-stream level in winter 



time as already observed above San Pedro Martir. Also a 
similar inversion in the wind speed gradient in the high part 
of the atmosphere above the jet-stream (at ~ 17 km a.s.l.) 
is visible above Mt. Graham as well as above San Pedro 
Martir in summer. Above Mt. Graham the median profile 
of the potential temperature shows that the typical 'slope 
change ' identifying the position of the tropopause visible at 
10 -11 km in winter time, is a little smoother and located 
at higher heights (~ 13 km a.s.l.) in summer. The level of 
the thermal stability plays obviously an important role to 
determine the strength of the optical turbulence triggered 
by the dynamic shear typical of the jet-stream. This fact 
might therefore play a role in the mechanism producing the 
shift towards higher heights of the jet-stream C% peak but 
is, in any case, coherent with what found above San Pedro 
Martir. The physical model proposed by Masciadri & Egner 
(2006) to explain the a effect seems therefore to be consis- 
tent with what is observed above Mt. Graham in this paper. 
We note that, the fact that the wind speed is visibly one of 
the main causes triggering the seasonal variation of the C% 
at the jet-stream level does not mean that for equivalent 
wind speeds above whatever astronomical site, one has to 
expect the same C% value. The absolute strength of the 
C% depends in a precise region of the atmosphere, indeed, 
on the thermodynamic stability of the atmosphere in the 
same region. 

6.2 HVR-GS: vertical distribution for h < 1 km 

The optical turbulence vertical distribution with high 
resolution (20-30 m) in the first kilometer is obtained 
with the method called HVR-GS presented in Egner & 
Masciadri (2007). The HVR-GR data-set is obtained taking 
the integral of the C% profiles retrieved from the AC 
frames and redistributing the energy in the first kilometer 
according to the detected triplets in the CC frames. Three 
different strategies can be used to study the turbulence 
spatial distribution in the boundary layer. The usefulness 
of each method depends on the application one intends to 
give to the analysis. We study: 

(A) the median of the C% (and/or J) profiles. 

(B) the average of the C% (and/or J) profiles. This 
method is very useful for comparisons of measurements and 
simulations obtained with atmospherical models. Moreover 
this operator has the advantage that the mean of Ji is equal 
to the J retrieved from the mean of the C% j (where Ji and 
C 2 N i refer to the each individual profile). It is not the case 
for the median. 

(C) the composite profiles as calculated in Section [5] This 
method is very useful for applications to Adaptive Optics. 

Following the strategy (A) we find that C% = for h 
> 25-30 m. This can be explained with the fact that each 
vector {C% profile) has many zeros. This is due to the fact 
that, during the monitoring, the system detects spikes and 
thin layers that changes position and duration with the time. 
In other words, the use of the median in these cases, can be 
misleading and it can provide very low values of the C% . 
Considering that the strategy (A) is not really useful to char- 
acterize these measurements in our case, the strategies (B) 
and (C) have been used. Figure [BJ-left shows the result ob- 
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Median CN2 Profile 
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C N 2 ( m - 2/3 ) C N 2 (m-^) 

Figure 4. Median profile obtained with the complete sample of 43 nights, the summer [April-June] and winter [October-March] 

time samples. Results are obtained with the standard GS technique. 



Potential temperature 



Median Wind Speed 




350 400 450 

Potential temperature (degrees) 

Figure 5. Left: Median potential temperature calculated using the ECMWF analyses. Right: Median wind speed calculated using the 
composite profile (GS below 2km, ECMWF analyses above 2 km - see text) in summer (thin style line) and winter (bold style line). 



tained following the strategy (B) i.e. the mean of the C 2 N 
profiles (Ah ~ 25-30 m) calculated from the sample 'HVR- 
GS' after normalization for the f 9 ; factor and it is therefore 
representative of 43 nights. The two profiles (with and with- 
out dome contribution) are shown in proximity of the ground 
(Figure[6]-right). To retrieve the typical scale height B of the 
exponential decay of the mean C% profile the measurements 



done below 125 m have been fitted with an exponential law 
(Eq.|T5|) as we already did in Stoesz et al. (2008): 

y = A ■ e ( -' l/s) (12) 

where A and B are free parameters. The calculation is 
obviously shown only in the case in which the dome contri- 
bution is subtracted. The fit gives A = 3.34-10" 15 and B = 
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37.4 m (Fig[6]-left). If we limit the analytical fit to the first 
30 meters, the scale height B = 28 m. The exponential decay 
is however an absolutely arbitrary analytical law. The im- 
portant issue if that these results definitely indicate that the 
HVR-GS technique is able to put in evidence that the turbu- 
lence decays above typical astronomical sites in stable night 
time conditions in a much sharper way than what has been 
predicted and quantified in the past. Indeed, the Hufnagel 
model in proximity of the surface (Roddier, 1981) in night 
conditions states that the turbulence scales as h -2 / 3 (Fig[6]). 
For the Hufnagel model (Fig[6]4eft) the C% decreases of one 
order of magnitude within 1 km while our results indicate 
that the C% decreases of one order of magnitude within ~ 
60-70 m (Fig[6]-right). It remains interesting the (800 m - 1 
km) range in which visibly a weak turbulence develops. One 
should expect a smoother connection between the boundary 
layer and the free atmosphere. We think that this is just an 
artifact effect due to the use of a quite different resolution 
below and above 1 km. From the point of view of the AO 
simulations this small gap should not cause any problems. 
It should be enough to implement a weak convolution to 
slightly smooth out the C% vertical profiles at the interface 
located at 1 km so to obtain a less abrupt connection of the 
C% profile above and below 1 km. 



We note that the presence of no signal (therefore C% = 
0) might potentially indicate that the turbulence strength 
is not equal to zero but simply weaker than the threshold 
C% =10~ 16 (associated to an equivalent J=2.5T0~ 15 ). In 
Masciadri et al. (2009) it has been calculated, in a post- 
processing phase, the most conservative case in which we 
assigned C% =10 -16 where there is no signal. After a more 
careful investigation we observed that the lattest distribu- 
tion is associated to a too large total J in the boundary layer 
and, for this reason, it can be discarded. 



/ h \ 

J C 2 N (h*)dh* I 



P(h) 



Figure [7] shows the results obtained following the strat- 
egy (C) that is the mean C% profiles associated to the J 
(or r ) related to the 20-30 %, 45-55 % and 70-80 % of the 
cumulative distribution. In Annex C the numerical values 
for the correspondent J values are reported. Curiously the 
first grid point near the ground of the 75% case distribution 
shows a weaker value with respect to the 25% and 50% cases. 
This is due to the fact that the third quartile dome seeing 
(e<2,75), that has been subtracted from the original value is 
particularly large (0.70 arcsec). The morphology of the C% 
in the first kilometer is very interesting showing several thin 
layers and a very weak turbulence between 800 m and 1 km. 
We note that the sample on which we calculate the aver- 
age in each slot [20-30]%, [45-55]% and [70-80]% is of a few 
hundreds of C% profiles therefore there are no doubts that 
this structure reproduces some real distrubution. In terms of 
morphology of the turbulence profile we find therefore that 
the higher the vertical resolution the thinner the size of the 
detected layers. This conclusion is perfectly coherent with 
the turbulence structure resolved by balloons equipped for 
the C% measurements (Azouit & Vernin, 2005). 



Figure [8] shows that the percentage of turbulence P(h) 
developed in the (0, h) range with h in the (0, 1 km) vertical 
slab. The function P(h) is defined as: 



JC 2 N (h*)dh* 



x 100 



(13) 



where the C N profile is that associated to the 45-55 
% case (Table IC1I and Fig. 0-centre). It is worth to note 
that the error bars for the HVR-GS technique is of the same 
order of half of the vertical resolution i.e. ± 12-15 m. This 
derives mainly by the definition of the zero point i.e. the 
ground that is characterized by an uncertaintly equivalent 
to half of the vertical resolution. 

From Fig[8] we retrieve that around 50% of the turbu- 
lence developed in the whole 20 km is concentrated below 
80±15m. This result is substantially different from the pre- 
liminary indications obtained in Egner & Masciadri (2007) 
and the turbulence seems to be much more concentrated 
near the ground. We also note that the C% morphology 
below 1 km decreases in a different way if we look at the 
sample 'WB' and the 'HVR-GS'. This is not a contradiction 
and it can be explained with the fact that the turbulence 
spatial distribution for the 'WB' and the 'HVR-GS' sam- 
ples is not necessarily the same because the first is based on 
the AC frames while the second is based on the CC frames. 
It is interesting to note that the higher the vertical resolu- 
tion, the sharper is the exponential decay of the morphologic 
turbulence structure. 

We note that, simultaneously to our study, some other 
authors (Chun et al. (2009)) recently investigated the turbu- 
lence structure near the ground at high vertical resolution. 
Even if the instruments employed were different these stud- 
ies present many similarities in the results. From a quali- 
tative point of view, we note that, also in that case, the 
turbulence appears well confined near the surface. From a 
quantitative point of view, things are more delicate. Look- 
ing at Table 3 (in that paper) it appears that their data- 
reduction is more similar to our method (C) than to the 
others methods. Figure [5] shows the J-profile retrieved from 
Fig. [5] - centre overlapped to the J-profile retrieved from 
Table 3 (Chun et al. (2009). While above Mt. Graham the 
turbulence vertical sampling is 25 m, above Mauna Kea the 
sampling increases from 15 m up to 80 m and it extends 
only up to 650 m. The turbulence vertical distribution ap- 
pears very similar, we note that above Mauna Kea a local 
minimum is present at ~ 45 m above the ground more or 
less in correspondence of the abrupt detection break due to 
the sensitivity threshold from LOLAS (Table 3 - Chun et 
al. (2009)). We highlight that the evident huge 'turbulent 
vacuum zone' between 560 m and 1 km (SLODAR) simply 
means that the turbulence is not measured in this vertical 
slab, not that turbulence is not present. 



7 CONCLUSIONS 

In this paper we present the results of a study aiming 
to characterize the optical turbulence at Mt. Graham. We 
present a general overview of the statistics (43 nights) of 
the Cjv profiles and all the main integrated astroclimatic 
parameters and their seasonal trends. The main conclusions 
we achieved are: 
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Figure 6. Left: Mean profile (A/i~20-30 m) calculated with the sample 'HVR-GS' and the dome seeing contribution included 

(black line). Hufnagel model (red line). Right: Zoom of the first 200 m. Near the ground it is put in evidence the different shape 
in case the dome contribution is included (black line) and excluded (green line). 
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Figure 7. Mean profiles calculated from the corresponding J (or ro) values related to the 20-30 %, 45-55 % and 70-80 % ranges of 
the J cumulative distribution. The dome contribution is excluded. In Annex C the numerical values for the correspondent J values are 
reported. 



(1) With a median seeing e= 0.95 arcsec (e= 0.72 arc- 
sec without dome contribution), isoplanatic angle 9q= 2.5 
arcsec and a wavefront coherence time to=4.8 msec, Mt. 
Graham confirms its good quality in terms of turbulence 
characteristics typical of the best astronomical sites in the 
world. All the integrated astroclimatic parameters (the see- 
ing, the isoplanatic angle, the wavefront coherence time and 
the equivalent wind speed) show a clear seasonal trend that 
indicates better turbulence conditions and weaker equiva- 
lent wind speed Vo in summer with respect to the winter. 

(2) The ground layer is characterized for the first time with 
a high resolution (200-250 m and 20-30 m). The turbu- 
lence exponentially decays above Mt. Graham with a much 
sharper profile than what has been supposed so far and ex- 
pressed with the Hufnagel model. Three different strategies 
of analysis aiming to investigate the morphology of the tur- 
bulence spatial distribution have been presented. We find 
that around 50% of the turbulence developed in the whole 
atmosphere is concentrated below 80±15 m from the ground 
and 60% of the turbulence in the first kilometer. This ev- 
idence together with the favorable large 9o observed above 
Mt. Graham (particularly in the spring/summer time) repre- 
sent extremely favorable conditions for astronomical obser- 
vations assisted by a LGS/GLAO system such as ARGOS. 

(3) We observe that the higher is the vertical resolution of 



the tool used to measured the turbulence vertical distribu- 
tion the sharper is the turbulence decreasing. 

(4) The percentage of turbulence developed below the pri- 
mary mirror of the LBT (i.e. ~ 35 m from the ground) is 
around 33%. However this estimate has to be considered 
with precaution because the uncertainty (2u~ 25-30 m) is 
of the same order of magnitude of the vertical resolution 
(AH ~ 25-30 m) and in the first hundred meters the turbu- 
lence decreases very sharply. 

(5) It appears evident that at Mt. Graham the turbulence 
decreases above the ground similarly to what observed above 
Mauna Kea in a more or less simultaneous study (Chun et 
al. 2009) performed with different instrumentation. 

(6) A composite distribution of the turbulence on the whole 
20 km is calculated to be used as input of AO simulations of 
the LBT Laser Guide Star system named ARGOS and the 
calculation of the gray zones for the near-infrared J, H and 
K band is done. The gray zone extends from a minimum 
of (200m - 378 m) in J band with the 'bad'-'bad' case up 
to a maximum of (200m - 3777 m) in K band in the case 
'good'-'good' case. 

(7) A clear C% seasonal variation trend has been observed 
in proximity of the ground and in the jet-stream regions. 
These measurements confirm the first evidence of the C% 
seasonal trend observed by Masciadri & Egner (2006) above 
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Figure 8. Percentage of turbulence developed between the ground and the height h with respect to the turbulence developed in the 
whole atmosphere (~ 20 km) as retrieved from the HVR-GS measurements and extended to the first kilometer. On the left hand is shown 
a zoom of the picture centred on the first hundreds of meters. 
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Figure 9. Red line: turbulence vertical distribution (J-profilc, 
50% case) above Mt. Graham as calculated in FigfTl-centre. Black 
line: turbulence vertical distribution (J-profile, 50% case) above 
Mauna Kea as calculated by Chun et al. (2009) in Table 6-centre. 
Mauna Kea measurements extend up to 650 m. 



other astronomical site. The physical model proposed by 
Masciadri & Egner (2006) that is able to explain the sea- 
sonal effect of the secondary peak of the C% called 'a ef- 
fect', is confirmed and refined. 

(8) For the first time we observed a seasonal trend of the 
dome seeing. This is certainly a topic that deserves a more 
careful investigation in the future. 

(9) We proved that the error in the normalization of the 
scintillation that has been recently put in evidence in the 
principle of the GS technique affects these measurements 
with an absolutely negligible quantity (0.04"). In other 
words, the median seeing retrieved from the GS (without 
correction) overestimates for ~ 0.04" the correct median see- 
ing. From a general point of view, all the GS measurements 
obtained with a pupil size D > 1.5 m and a binary separa- 
tion 9 < 8 arcsec, are affected by this error of less than a 
few hundredths of arcseconds. 



SCINTILLATION FOR THE STANDARD 
GENERALIZED SCIDAR 

It has been highlighted (Johnston et al. 2002) an impreci- 
sion in the normalization of the autocovariance of the scin- 
tillation maps obtained with a GS. A generalization of this 
problem, extended to altitudes h > 0, has been presented 
more recently by Avila & Cuevas (2009). Results of these 
studies say that, to obtain exact results for the C% at an 
height h from the ground, one has to multiply the C% (h) 
retrieved from the GS by a factor l/(l+e(/i)) where e is the 
relative error between the exact and erroneous autocovari- 
ance of the scintillation maps obtained with the GS. e(h) 
depends on a set of geometrical parameters related to the 
optical set up and the observed binaries, more precisely the 
pupil of the telescope D, the height d at which the detection 
plane is conjugated below the ground, the ratio between the 
stellar magnitude of the binary, the angular separation 9 and 
the central obscuration of the pupil e expressed as a fraction 
of the primary mirror. The authors introduced an approxi- 
mated solution for the correction valid for a simple circular 
telescope pupil and the exact solution valid for a telescope 
pupil shape formed by a primary and a secondary mirror. 
However we calculated that, for our sample, the correction 
of all the C% profiles has an absolutely negligible impact on 
the statistical analysis presented in this paper (in both cases 
i.e. approximated and exact solution) and the median e cal- 
culated on the whole 43 nights sample changes respectively 
for just 0.03" and 0.04". This means that approximated and 
exact solution are very similar in this case. 

To prove that, we first note that in the LBT-GS, the 
scintillation maps and the corresponding autocovariance cal- 
culated in real time are sampled on the same number of pix- 
els square. If the pupil of the telescope takes the whole field 
on the CCD, under the condition we have just described, it 
follows that hmax = D/(2-f?) (and not D/f? as indicated by 
Avila & Cuveas (2009)). It is known, indeed, that a layer at 
a height h produces a couple of identical scintillation maps 
separated by a distance d=9-h at which are associated the 
later peaks of the auto-correlation located at a distance d 
from the centre of the frame. However, the most external 
lateral autocorrelation peaks visible on the CCD (that cor- 
respond to the highest detectable layer) are necessarily lo- 
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Figure A2. Median profiles related to 43 nights (black line) 
and corrected by the normalization factor (red line). 



cated at a distance d—D/2 (i.e. half of the CCD size) from 
the centre of the frame if the pupil of the telescope covers the 
whole CCD and if the support of the auto-correlation is not 
changed. It follows therefore that h max — D/(2-#). The later 
peaks associated to turbulent layers located at h > D/(2-#) 
fall outside the CCD. This means that the vertical range 
monitored by the LBT-GS is basically half than what is in- 
dicated by Avila & Cuevas (2009). The important condition 
to be respected is that h max includes the whole atmosphere 
and this is the case: the smallest h max is 24 km in the case 
of the widest binary (7.6 arcsec). We are therefore in a cor- 
rect configuration to scan the whole atmosphere. Besides we 
also note that, as shown in Masciadri et al. 2002-Fig.l2 for 
example, the statistical noise increases for h > 20 km (at 
least for the geometry commonly used by the standard GS 
and D > 1.5 m) and the GS measurements are no more re- 
liable in that region. Looking at Fig. 4 and Fig. 5 (Avila & 
Cuevas, 2009) it is possible to appreciate that, for D > 1.5 
m and e = 0.2 (as it is our case) and for b = a/(l+ a) 2 , a 
= io-° 4Am , the value of e(h) in the exact solution (Sec. 3.2 
of that paper) and the approximated one (Sec. 3.1 of that 
paper) is very small in both cases in the (0, h ma x) range. 

Figure IA1I shows the relative error e as a function of h 
for all the different binaries used in this study (Table [3J) and 
for a fixed value of D = 1.83 m, e — 0.2, d = 3.5 km for 
the exact solution (full line) and approximated one (dashed 
one). We observe that the error e, in both exact and approx- 
imated cases, is typically smaller than 0.1 in the 20 km range 
and it reaches the value of 0.15 only with the widest binary 
7 Ari (6 = 7.6"). The error produced on the C% recon- 
structed by the normalization of the autocovariance of the 
scintillation maps is absolutely negligible in our case as can 
be seen in Fig l A2I producing a difference in the seeing of the 
order of 0.03" and 0.04" respectively with the approximated 
and exact solution. 



APPENDIX B: NORMALIZATION OF 
SCINTILLATION FOR THE HVR-GS 

The error in the normalization put in evidence by Johnston 
et al. (2002) and Avila & Cuevas (2009), in spite of the 
wider binary employed, affects the HVR-GS measurements 



in the same negligible way in which the standard GS tech- 
nique is affected. The reason is that this technique, thanks 
to the procedure called 'normalization' (see Sec. [SJ), uses 
the AC and the CC frames just to retrieve the vertical tur- 
bulence distribution in the atmosphere (i.e. the shape). The 
strength of the turbulence is corrected by the factor f 9 ; that, 
physically speaking, is equivalent to taking the turbulence 
measured with the standard GS and redistributing it follow- 
ing the shape of the profile reconstructed by the HVR-GS 
technique. The HVR-GS technique retrieves the C% pro- 
files in the first kilometer. Let's assume to have an error A 
= l/(l+e) for the C% as indicated by Johnston et al. (2002) 
and Avila & Cuevas (2009). We should therefore expect that 
the corrected turbulence strength is Jo * = A -Jo for each 
vertical grid point of ~ 25 m. However, when we multiply 
for the f 9 ; factor fEq II 1 [I . the factor A disappears. Indeed f 9 ; 
can also be written as Jo /(A-Jo ). We conclude that, what- 
ever is the value of A, the error potentially introduced by 
the wide binary is not taken into account in the final value 
of the C% retrieved with the HVR-GS technique. 



APPENDIX C: VERTICAL DISTRIBUTION 
FOR THE HVR-GS: J VALUES 

In Table [Cll are reported the values of J = C% -dh obtained 
with the average of J included in the 20-30%, 45-55% and 
70-80% ranges. 
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Table CI. Composite profiles for h < 1 km after the 'normalization' for the f g ; factor and the dome contribution included. These 
composite profiles are statistically representative for 43 nights. In the first column the height, in the second, third and fourth columns 
the J values for the 'good', 'typical' and 'bad caes. In the last line are reported the J values without the dome contribution, (median 
values: £<i,25 = 0.35 arcsec, £<j,50 = 0.52 arcsec, £<j,75 = 0.70 arcsec). The first grid point (h = 0) includes Cf^ values between -12.5 m 
and +12.5 m. 
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